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ABSTRACT 

We present an extension of the Kompaneets equation which allows relativistic 
effects to be included to any desired order. Using this, we are able to obtain 
simple analytic forms for the spectral changes due to the Sunyaev-Zel'dovich 
effect in hot clusters, correct to first and second order in the expansion parameter 
Oe = ksTe/mc 2 . These analytic forms agree with previous numerical calculations 
of the effect based upon the multiple scattering formalism, and are expected to 
be very accurate over all regions of the CMB spectrum for k B T e up to ~ lOkeV. 
Our results confirm previous conclusions that the result of including relativistic 
corrections in the Sunyaev-Zel'dovich effect is a small reduction in the amplitude 
of the effect over the majority of the spectrum: specifically we find AT/T = 
—2y(l — 17/10 9 e + 123/40 # e 2 ) (correct to second-order) in the Rayleigh- Jeans 
region, where y is the usual Comptonization parameter. For a typical cluster 
temperature of 8keV, this amounts to a correction downwards to the value of 
the Hubble constant derived using combined X-ray and Rayleigh- Jeans Sunyaev- 
Zel'dovich information by about 5 percent. 



Subject headings: cosmology: cosmic microwave background — galaxies: clusters: 
general — radiative transfer — scattering 



1. Introduction 

Non-relativistic treatments of the Sunyaev-Zel'dovich effect usually employ the Kom- 
paneets equation |(Kompaneets 1957)| to determine the distortion of the Cosmic Microwave 
Background (CMB) spectrum. The Kompaneets equation does not, however, include rel- 
ativistic effects, which may be important for hot clusters where ksT e ^ lOkeV. For this 
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reason, and because of the low optical depth of typical clusters, relativistic treatments of 
the Sunyaev-Zel'dovich effect usually employ a multiple scattering description of the Comp- 
tonization process |( Wright 1 979)1 KFabbri l98l)| , |(Taylor fc Wright 1989)| , |(Loeb, McKee^j 



Lahav 1991) , (Rcphacli 1995) . Including relativistic effects in this procedure gives a compli- 



cated expression for the spectral distortion, which is best handled by numerical techniques 
(see, for example, |(Rephaeli 1995)|) . 

In this paper, we show how the Kompaneets equation may be extended to include 
relativistic effects in a self-consistent manner, allowing the Sunyaev-Zel'dovich effect in hot 
clusters to be described on the basis of a Kompaneets type equation. The extension of the 
Kompaneets equation can be carried out to arbitrary orders in relativistic effects, although 
we shall only consider the lowest order corrections here. The resulting equation conserves 
photons at every order, and we demonstrate that it is consistent with earlier calculations of 
the energy transfer rate between the plasma and a Planck distribution of photons (Woodward 



1970 



We then consider the application of the generalised Kompaneets equation to the calcula- 
tion of the Sunyaev-Zel'dovich effect in hot clusters. Simple analytic forms are given for the 
spectral distortions in the limit of small optical depth, including relativistic effects to second- 
order. These are in excellent agreement with Rephaeli's (1995) numerical calculations, which 
were based on the multiple scattering approach (truncated at one scattering). This lends 
further support to Fabbri's observation that the Boltzmann equation can be applied to de- 
scribe the Sunyaev-Zel'dovich effect in optically thin clusters |(Fabbri 1981 )| , despite claims 
to the contrary |( Wright 1979) . 



The relativistic corrections to the Sunyaev-Zel'dovich effect are of importance in the 
calculation of the Hubble constant H by the Sunyaev-Zel'dovich route in hot clusters (see, 
for example, |(Lasenby fc Jones 1997)| for a recent review of the Sunyaev-Zel'dovich route of 



determining Hq, and |( Saunders 1996 )| for recent observations). In the Rayleigh- Jeans region, 
we find that relativistic effects lead to a small decrease in the Sunyaev-Zel'dovich effect, 
and hence a small reduction in the hitherto determined values of H Q , in agreement with the 



conclusions in |( Rcphacli fc Yankovitch 1997) 



We employ natural units, c = Ti = 1, except in the discussion of the Sunyaev-Zel'dovich 
effect in Section f|. 
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2. Extending the Kompaneets Equation 



In this paper we shall not consider effects due to the peculiar motion of the cluster 
(such effects give rise to a kinetic correction to the Sunyaev-Zel'dovich effect). For a comoving 
cluster, the CMB photon distribution function is isotropic and may be denoted n(u), where u 
is the photon energy. The electrons are assumed to be in thermal equilibrium at temperature 
T e , and are described by an isotropic distribution function f(E), where E is the photon 
energy. The Boltzmann equation describing the evolution of n(u) may be written as ( Buchlcr 
|& Yeuh 1976) 



dn{ui) 
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where W is the invariant transition amplitude for Compton scattering of a photon of 4- 
momentum k^ by an electron (of charge e and mass m) with 4- momentum p^, to a photon 
momentum k'^ and an electron momentum p ,fl |(Berestetskii, Lifshitz, fc Pitaevskii 1982)] : 
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with k = —2p 1 k l j l and k' = 2p^k'^. In equation ( |2- 1[ ) , we have assumed that electron degen- 
eracy effects may be ignored. 

The electrons are described by a relativistic Fermi distribution. Since we are ignoring 
degeneracy effects, we have 
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Substituting this form for f(E) into equation ( p-l|) , and expanding the term in brackets in 
the integrand in powers of Ax, where 



to 



X 



Ax 



gives a Fokker-Planck expansion 
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where 

In = ~J ^- 3 d 3 p'd 3 k'Wf(E)(Ax)\ (2-8) 

TV. J ( Z7T ) 



which does not depend on n(w). 



The calculation of the I n may be performed by expanding the integrand in powers of 
p/m and u/m. The factor f(E) is handled by the expansion 

f(E) « ~ "OAn^-u^ + i^ M 2 + i6) e V( u - 4) + • ■ ■), (2-9) 

where 

2mkBT e ^ ^ 

e e = (2-ii) 

m 

and the chemical potential /x may be eliminated by introducing the electron number density 
N e , which evaluates to 

N e = !^^T(3/2)(k B T e )le^- m V kBT e (l + ^+ ^ + 0(6!)) . (2-12) 



7T 2 v ' /v ' V 8 128 

Note that equation ( |2-12| ) is an asymptotic expansion about 9 e = of the result 

7T 



iV e = ^-^ 2 (l/0 e )e^/ fc ^, (2-13) 



where ^(a?) is a modified Bessel function, which suggests that our series expansions of the 
J n will only be asymptotic series. The calculations of the I n are ideally suited to symbolic 
computer algebra packages (we use Maple). Expressing the results in terms of 9 e , x, and the 
Thomson cross section ay, we find 

I x = \o T NAx ((4 -x) + 6> e (10 - f x + fx 2 )) + 0(6> e 3 ) 

h = \a T N e 6 e x 2 (l + 6 e (%-fx + l Q x 2 ))+0{6*) 

h = ±a T N e d e x 3 (l(G-x)d e )+0(63) 

h = \a T N e e e x\7e e ) + 0(6*). (2-14) 

For n > 4, I n is third-order or higher in 9 e . For CMB photons passing through a cluster at 
redshift z, we have x ~ 6.2 x 10 _4 (1 + z)//csT e , where x is the average of x over a Planck 
distribution and fc_eT e expressed in eV. The electron temperature is typically ^ lOkeV, so 
that i C 1. 
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Substituting these expansions for the I n into the series ( |2-7| ), we develop an expansion of 
dn/dt in 9 e . In this paper, we shall only be concerned with the lowest order relativistic cor- 
rections, and start by retaining all terms up to 0(# e 2 ). To include all such terms consistently, 
it is necessary to retain only the first four terms in the series fl2-7| ). A lengthy calculation 
gives the result 



dt 

where the current j(x) is given by 
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The zero-order term in equation ( 2-16 ) is just that term which usually appears in the Kompa- 



neets equation |(Kompaneets 1957)1 The 0(6 e ) term is the lowest-order relativistic correction 
to the current. The form of equation ( |2-15|) ensures conservation of the total number of pho- 
tons, which is true for each order in 9 e . A similar expression has been derived independently 
by Stebbins (1997) using non-covariant methods, although he only considers the lowest-order 
terms in u/m. 

This derivation makes it clear that trying to include relativistic terms in the framework of 
the usual Kompaneets equation |(Kompaneets 1957]! , although it would seem natural to retain 
only terms in I\ and I2, would be incorrect, since it would not obey photon conservation. To 
be specific, if only I\ and I2 are retained, then photon conservation requires that the current 
be given by 

. , s 1 2 f dn 
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and consistency with the expansion of the Boltzmann equation 
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Calculating I 2 to 0(6^) in the limit of small x, gives the current j'(x) (in the limit of small 



x) where 



j (x) = x 




(2-19) 



which amounts to a relativistic correction to the cross section appearing in the usual Kompa- 
neets equation. This result is incorrect, since the approach is not self-consistent to the order 
that the answer is quoted. To see this, one need only calculate I\ using equation ( |2-18| ) and 
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the expression for I 2 from equation ( 2-14 ). The expression for 1% obtained by this procedure 
only agrees with the direct calculation equation ( |2-14| ) to 0(6> e ), not to O(0 e 2 ). One obtains 
a better approximation to the true current (eq. ( |2-16|) ) by calculating I\ directly and then 
using equation (|2-18;) to calculate an effective value of which may then be used to deduce 



the current. This procedure reproduces the first term in the 0(6^) correction in j(x), in the 
limit of small x, and is similar to that employed by Fabbri (1981) in obtaining his equation 
(11). However, it is only the full expression equation (|2-16|) which is consistent with the 
original Boltzmann equation (|2-1|) to 0(6^) 



3. Rate of Energy Transfer 

As a check on the consistency of equation ( |2-l(j ) with existing results in the literature, 
we calculate the energy transfer rate between the electrons and a Planck distribution of 
photons at temperature T r . 

Multiplying the continuity equation ( [2-1 5| ) by x 3 and integrating, we find 

Q rco roc 

— / x 3 n(x) dx = —aTN e 6 e \ x 2 j(x)dx. (3-1) 
ot Jo Jo 

The left-hand side of equation ( |3-1[ ) is proportional to the rate at which the photons are 
gaining energy per unit volume, denoted dE /dt. Substituting a Planck distribution for n(uj) 
in equation (|2-16| ) and integrating gives the result 

^ = 4E r N e a T (6 e - 6 r ) (l + h e - 2l||a + 0(# 2 )) , (3-2) 

where 

Br = (3-3) 
m 

and ((x) is the Riemann Zeta function. 

This result may be compared with a direct evaluation of the energy transfer, obtained 
by multiplying the transition rate W by the energy transfer uj' — u, and integrating over all 
collisions: 

^T = 4 / H + n{u'W - u). (3-4) 

The integral may be evaluated by a consistent expansion of the integrand. A lengthy calcu- 
lation gives the result (|3-2|) , which also agrees with the non-covariant calculation of Wood- 
ward (1970) (who gives higher-order corrections also). Note that if we had considered only 
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Ii and then using j'(x) given by equation ( |2-19|) , we would have obtained (in the limit 
of small 9 r ) 

d,E / 47 \ 

= AE r N e o T {6 e - %) (l + -^9 e + 0{6 2 e ) J , (3-5) 

which grossly overestimates the importance of the relativistic corrections. 



4. The Sunyaev-Zel'dovich Effect 



In this section we apply the generalised Kompaneets equation (to first-order in relativis- 
tic corrections) to the calculation of the Sunyaev-Zel'dovich effect in optically thin clusters. 
We consider higher-order effects in the next section. 

Following the standard assumptions, we assume that the optical depth is sufficiently 
small that the spectral distortions are small. In this limit, we may solve equation ( |2-15| ) 
iteratively. The lowest order solution is obtained by substituting the initial photon distribu- 
tion no(x) into the current (eq. ( |2-16 )). The integral over time is then trivial, and may be 
replaced by an integral along the line of sight through the cluster, giving 

M*) = (ftW) , (4-1) 

where j(x) is evaluated with n (x), and 



y = a T N e e dl 



(4-2) 

where the integral is taken along the line of sight through the cluster. 

For the CMB we take the initial (undistorted) photon distribution to be Planckian with 
temperature T : 

(4-3) 



n (x 



- 1 



where a = T e /T is the (large) ratio of electron temperature to the CMB temperature. 
Evaluating equation (|4-1|) in the limit of large a, we find the following fractional distortion: 
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correct to first-order in relativistic effects, where 
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9 e AI/y 




Fig. 1. — The intensity change 6 e Al/y (in units of 2(A;sT ) 3 /(he) 2 ) plotted against X for 
three values of ksT e (in keV). The solid curves are calculated using the first-order cor- 
rection to the Kompaneets equation, while the dashed lines are calculated from the usual 
Kompaneets expression. 



The first two terms in square brackets in equation ( |4-4| ) give the usual non-relativistic 
Sunyaev-Zel'dovich expression, while the terms proportional to 8 e are the lowest order rela- 
tivistic correction. Equation ( |4-4p agrees with the result in Stebbins (1997). In the Rayleigh- 
Jeans limit (small X), we find 

A "M~ a »('l-iZft + 0«f)V (4-6) 



n(X) * V 10 



In Figure [I] we plot the change in spectral intensity 6 e AI /y as a function of X, where 

X 3 An 
e x - 1 n 



Also plotted in Figure [l] are the non-relativistic predictions made with the standard Kompa- 
neets equation. The curves in Figure p] are for ksT e = 5, 10, and 15keV, which are the same 
as the parameters used by Rephaeli (1995) in his Figure 1. His calculations, which were 
based on the multiple scattering formalism |( Wright 1979)| and required a numerical analysis, 
give results in excellent agreement with ours, which only require the use of the simple expres- 
sion (|2-8|) . This suggests that there is no problem in principle with applying the Boltzmann 
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Fig. 2. — The fractional change AT R j/To in the Rayleigh- Jeans brightness temperature 
plotted against X for three values of k B T e (in keV). The solid curves are calculated using 
the first-order correction to the Kompaneets equation, while the dashed lines are calculated 
from the usual Kompaneets expression. 

equation to the problem of Comptonization in clusters even though the optical depth may 
be very small. Similar conclusions were reached by Fabbri (1981), but his demonstration 
was restricted to low temperature clusters where relativistic effects are not important. 

It is clear from Figure |I] that for X <^ 8, the relativistic corrections lead to a reduction 
in the magnitude of the intensity change, compared to the non-relativistic prediction. This 
in turn leads to a reduction in the inferred value of the Hubble constant determined by the 
Sunyaev-Zel'dovich route. 

In Figure ^| we plot the fractional change in the Rayleigh- Jeans brightness temperature 



for the same 6 e as in Figure [T| The relativistic corrections to the change in the Rayleigh- Jeans 
brightness temperature are significant even at low frequency, unlike the corrections to the 
intensity, where relativistic corrections are small in the Rayleigh- Jeans part of the spectrum. 



AT RJ /T (divided by y/9 e ), where 
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5. Higher-order Effects 

We have found that for k B T e ^ lOkeV the second-order relativistic effects make a 
significant contribution to the spectral distortion, while third-order effects are only significant 
for k B T e £ 15 keV. 

These calculations require a straightforward extension of the method of Section [| to 
include terms at 0(6f) (for second-order relativistic effects). For the calculation to 0(# 3 ), it 
is necessary to retain the first six terms of the series (|2-7| ), and to calculate Ii through I 6 to 
0(6* e 3 ). The first iteration of equation ( 2-15 ) for T e ^> To gives the following next order (in 
6 e ) correction to An/n: 
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In the Rayleigh- Jeans limit, we find 

An(X) ^ 17„ , 123 o 



+ m ., A 4/i^ (-528 + 187Xcoth(|X))1. (5-1) 



n(X) -- 2 H 1 -T0 ft+ 4iT 9 « +O ^J- < 5 " 2 ) 

In Figure || we compare the spectrum of AJ calculated with equation ( [4-4j ) to the 
spectrum with the correction ( |5-1J) included, for ksT e = 5, 10 and 15keV (9 e ~ 0.01, 0.02 and 
0.03 respectively). In each case, the second-order relativistic effects are not significant in the 
Rayleigh- Jeans part of the spectrum. This is to be expected from inspection of equation ( |5-2| ) , 
where the # e 2 term is clearly insignificant for the values of 9 e considered. For ksT e = 5 keV, the 
second-order effects are insignificant over the entire spectrum. However, for k B T e ^ lOkeV, 
the second-order effects make a significant contribution to the relativistic correction to the 
Kompaneets based prediction outside the Rayleigh- Jeans region. We have verified that the 
third-order corrections are negligible over the entire spectrum for k B T e ~ lOkeV. This is 
confirmed by a comparison of the curves in Figure |3] with the points which are the results of 
a direct Monte-Carlo evaluation of the Boltzmann collision integral with n(u) given by the 
Planck distribution no{x) |( Gull fc Garrett 1997 )| . The second-order effects should be included 



in the analysis of high frequency data for hot clusters. The magnitude of the second-order 
correction to the Sunyaev-Zel'dovich result for the rather mild values of 9 e considered here, 
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Fig. 3. — The intensity change d e AI/y (in units of 2(A;#T ) 3 /(he) 2 ) plotted against X for 
three values of ksT e (in keVj. The solid curves are calculated using the second-order correc- 
tion to the Kompaneets equation, while the dashed lines are calculated from the hrst-order 
correction. The points are the results of a Monte-Carlo evaluation of the Boltzmann collision 
integral. 

is symptomatic of the asymptotic nature of the series expansion of dn/dt in 9 e . However, for 
the majority of clusters considered in Sunyaev-Zel'dovich analyses, the inclusion of the first 
two relativistic corrections should be sufficient, particularly for experiments working in the 
Rayleigh- Jeans region of the spectrum. 



5.1. The Crossover Frequency 

The accurate determination of the crossover frequency X (where the thermal com- 
ponent of the spectral distortion vanishes) is essential for reliable subtraction of the kine- 
matic contribution to the Sunyaev-Zel'dovich effect |(Kephaeli 1995)| . In Figure [| we plot 
the crossover frequency as a function of ksT e , with the first three relativistic corrections 
included. For fcsT e <^ 20 keV we find that X is well approximated by the linear relation 

X ~ 3.83(1 + 1.136> e ). (5-3) 



5 ' ' ' ' 10 ' ' ' ' 15 ' ' ' ' 20 

k B T e /keV 



Fig. 4. — The crossover frequency X plotted against k B T e . The solid line is calculated with 
the inclusion of third-order corrections to the Kompaneets equation. The upper dotted line 
is a linear fit to the solid line with X Q = 3.83(1 + 1.13# e ), while the lower dotted line is the 
linear fit given in \(Rephaeli 1995} : X = 3.83(1 + e ). 



For comparison, Rephaeli (1995), found X Q to be approximated by X ~ 3.83(1 + 6 e ) in the 
interval k B T e = 1-50 keV, while Fabbri (1981) found X ~ 3.83(1 + 1. 10 e ) for k B T e £ 150 keV. 
It is clear that our calculation favours Fabbri's expression. For k B T e ^ 20keV, X calculated 
with the first three relativistic corrections departs from the linear prediction ( |5-3| ). However, 
we do not regard this as indicative of a breakdown of the linear approximation, since it is 
clear from Figure |3| that the inclusion of higher-order terms may have a significant effect on 
the value of the crossover frequency. 



6. Conclusion 



We have shown how the Kompaneets equation may be generalised to include relativistic 
effects in a self-consistent manner. The resulting equation guarantees photon conservation 
at each order and is in agreement with direct calculations of the energy transfer rate between 
the plasma and a Planckian distribution of photons. 

We have applied this formalism to the calculation of the Sunyaev-Zel'dovich effect in 
optically thin clusters. We presented simple analytic expressions for the first two relativistic 
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corrections to the usual Kompaneets based expression for the spectral distortion, which 
are in excellent agreement with the numerical calculations of Rephaeli (1995) for electron 
temperatures <^ 10 keV. This provides further evidence that the low optical depth of clusters 
does not forbid the application of the Boltzmann equation to the calculation of the Sunyaev- 
Zel'dovich effect |(Fabbri 198 1 )| . The asymptotic nature of the expansion of dn/dt in 8 e 
requires the inclusion of higher-order corrections to calculate the effect in hotter clusters in 
the Wien region of the spectrum. While the calculation of higher-order corrections is not 
problematic, the bad convergence properties of the series means that ultimately one must 
resort to a numerical calculation of the collision integral |( Gorman 1970 )| (or employ the 
multiple scattering formalism |( Wright 1979)|) to calculate the effect in very hot clusters. 

Our calculations support the conclusions reached in |(Kephaeli fc Yankovitch 1997)| , that 
including relativistic effects leads to a small decrease in the value of the Hubble constant, 
inferred from combined X-ray and Sunyaev-Zel'dovich information. For a cluster temperature 
of ~ 8keV, the reduction in Hq due to relativistic effects is ~ 5 percent for measurements 
made in the Ray leigh- Jeans region. 
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